A New Class of Elliptic Finite-Gap 
Densities of the Polar Operator 
and Stationary Solutions of the Harry Dym equation 

L. A. Dmitrieva, D. A. Pyatkin 
Department of Mathematical 
& Computational Physics, St. -Petersburg State University, 
198504 St. -Petersburg, Russia 
e-mail: mila@JK1454.spb.edu 
danila@DP2627.spb.edu 

Abstract 

A new family of one- and two-gap elliptic densities of the polar operator has been 
constructed by modifying the so-called "higher time approach" to constructing finite-gap 
solutions of the Harry Dym equation. Auto-Backlund transformations of the obtained 
stationary solutions have been constructed and their properties have been studied. 

1 Introduction 

The present work is devoted to constructing elliptic one- and two-gap densities of the 
polar operator 

i— hw^, (ii) 

with periodic function r(x): r[x + T) = r[x) on the real axis x G R. Function r(x) 
is referred to as density of the polar operator. It is easy to see that the operator L is 
intimately connected with the string equation 

^ x + \p{x)^ = Q (12) 

provided that p(x) = r~ 2 (x) where p(x) is the density of the string and ip the amplitude of 
it's oscillation. We construct and study functional properties of some finite-gap densities 
of the operator (|Tl|) , namely 1-gap and 2-gap elliptic densitie^]) of the operator L. 

Despite classical character of the spectral problem (|l2|) for the operator L and a great 
number of papers devoted to the string equation (see works [||, [§], ||, j|, || and reference 
therein to name a just few), hitherto there weren't any explicit analytical formulae for 
the finite-gap densities r(x) with the exception of expression for smooth 1-gap density 
constructed explicitly in Here by explicit representation we assume representation of 
the form 

r(x) = R(y), 

where y = y(x) is solution of the functional equation 

x = F(y), 

L ) we call a density elliptic if it can be expressed in terms of elliptic functions. 
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and the functions R and F are given explicitly. 

The absence of analytical formulae for the finite-gap densities doesn't seem strange 
despite the relation between the string equation (12) and the one-dimensional Shrodinger 
equation 

-n; y + u ( y )n = xn (13) 

for which the class of finite-gap potentials is well known ||]. This relation has the form 
u(y) = y 2 rl(x) - 74 r(x) r xx (x) = (ft- 1 / 2 ) R 1/2 

\ / yy 



y - 1 7(kj 

if)(x,X) = Q(y,\)r 1/2 (x). 

So if the density r(x) is known then the potential u(y) can be reconstructed easily. How- 
ever the above relations can not be used for effective constructing the densities of the 
polar operator directly. 

At the same time the fact of finite-gap densities' existence was proved by M. G. Krein 
more than 40 years ago. In the paper Q a theorem was stated (without proof) asserting 
that there exist 2 N ~ l iV-gap densities of the string equation provided that the left edge 
of the operator's L zonal spectrum is situated at the point A = 0. The last condition 
in the Krein's theorem is essential)^), since in case of it's violation generally speaking 
the number of A^-gap densities is greater than 2 N ~ 1 . In the present paper we illustrate 
this effect by studying arrangement of the operator's spectrum in regard to the point 
A = and in case of 1- and 2-gap elliptic densities show that the set of all 1- and 2-gap 
densities is considerably wider and contains also some singular densities which were out 
of consideration in |jj . 

In order to construct smooth and singular 1- and 2-gap densities of the operator L we 
make use of the relation of linear spectral problem for the polar operator L with solutions 
of nonlinear evolutionary integrable Harry Dym equation [|To|]: 

r t + r 3 r xxx = 0, r = r(x, t), (14) 

and it's higher analogues.^) This relation is based on the fact that || holding the time 
variable t in the solution of the higher Harry Dym equation the resulting function turns 
out to be finite-gap density of the string equation. Because of this and the purpose of our 
paper (as we interested mostly in densities of operator L) we will not distinguish between 
densities of the string equation and solutions of Harry Dym equation. 

For the first time the finite-gap solutions of the Harry Dym hierarchy were constructed 
in H using higher times method developed in |7|. In |j] the finite-gap solutions were 
constructed in terms of multi-dimensional Riemann 0-functions and needed considerable 
effectivization. In the present paper we modify the scheme for construction finite-gap 
solutions of the Harry Dym equation and reduce the O-functions form in case of 2-gap 
densities to elliptic functions. This reduction becomes possible after the Riemann surface 
of the polar operator's Bloch functions is chosen to possess symmetry of the special sort. 



Along with higher times technique we use auto-Backlund transformation [11] for higher 
equations of Harry Dym hierarchy. We show that auto-Backlund transformation allows 
one to construct new 1-gap densities of the polar operator. We also construct auto- 
Backlund transformations in 2-gap case and study their properties leaving the property 
of spectrum preservation as hypothesis. 



2 - ) note that location of the Schrodinger zonal spectrum's edge in the point A = is not essential 



3) 



hierarchy of the Harry Dym equations is described later 



2 



2 Modified scheme for construction of finite-gap 
solutions of Harry Dym equation. 

The scheme for construction of iV-gap (N < oo) solutions of the Harry Dym equation 
proposed in || is the following. The iV-gap solution of the eq. (|TJ|) r(x,t\) is assumed 
to depend on N — 1 additional parameters — higher times t2, ■ ■ ■ , ijv in such a way that 
solution r(x,t) = r(x,ti, . . . ,tjy) with respect to t m , 1 < m < N, obeys the m-th higher 
Harry Dym equation: 

n m + r 3 {dlrlr) m r~ 3 r x = 0, m = l,...,N. (25) 

Here the operator I is defined by equations I d x r~ l = r _1 — 1 and I d x ip(x,t) = ip(x,t), 
if ip is any differential polynomial of r. 

The iV-gap solutions of the system (|25|) were obtained p| in the form: 



r(x, t) = R(y, t) = j±- ^ (26) 

l + (a,f(y,t)) 

where y = x + e(x,t). The phase function e(x,t) is related with solution r(x,t) by 
equation e x = 1/r — 1 and finally should be determined by solving functional equation 



y 

e(x,t) = E(y,t) = J(l-R(y',t))dy'. 



o 

In (^6|) by < • > we denote the inner product in H N . The constant vector a and the 
function f(y,t) have been described explicitly in Q in terms of iV-dimensional Riemann 
B-function and the periods of holomorphic differentials on the hyperelliptic curve 

2N 

A 2 = II( A -^)> E i^ E j (27) 

i=0 



with real numbers Ei coinciding with edges of spectral gaps of the polar operator (11) 
associated with Harry Dym equation. 

Remind that explicit expressions for f m (y, t) were obtained in || basing on the relation 

fm = £t m (x,t), 



and on the link |TT[ between Harry Dym and Korteweg-de Vries (KdV) hierarchies. 
Namely let u(jj,t\, . . . , tjv) be a solution of the first N equations of the KdV hierarchy, 
then 

u tm = 1 hd y e tm+1 , m = l,...,N. (28) 

In particular u = — 1 /2^t 1 - 

In the present paper we modify the scheme of construction of finite-gap solutions of 
the Harry Dym equation using Hamiltonian structure of the KdV hierarchy. Namely the 
system of the first N equations of the hierarchy 

U t m = -L m u y , m = l,...,N, 

where 

y 

L = dy — Au — 2u y J dy ■ 
can be expressed in Hamiltonian form: 



dy \2 5u 
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Here the functionals C m [u] represent the countable set of conservation laws of KdV equa- 
tion: 

oo 

C m = (-l) m (2m + I)" 1 J dy L m {y u y (y, t) + 2u(y, t)}, (210) 



and the symbol 5/5u in (29) denotes variational derivative: 



^ _ V 1 ( d \ 3 d (j) — ® 3u 



In case of periodic function u one should regard integrals in (|210| ) as integrals over pe- 
riod of u. But for us the most important fact is that from (f29|), (|210| ) it is readily seen that 
variational derivatives of conservation laws C m [u] appear to be differential polynomials of 
u(y,t). Comparing (e|) and ( ^9|) we obtain 



( _ ir «CpM + ^ (211) 

du 

where (3 m are constants which in periodic case can't be deduced from relations (^) and 
fl29|) . Computation of j3 m can be carried out in each particular case (1-gap, 2-gap, etc.). 

Thus the meaning of the scheme of solutions construction modification proposed in the 
present paper is in representation of f m (y,t) = Et m in the form of differential polynomial 
of u(y,t) by virtue of (|211|) . Then, knowing the solution of KdV equation possessing 
desired spectral properties we can compute all the necessary f m (y,t) and putting them 



in (26) obtain solution of the Harry Dym equation possessing the same spectrum. 



In particular from ( |211| ) and ( |210| ) it follows that 

f 1 {y,t) = e tl = -2u + l3 l (212) 
f 2 (y, t) = e t2 = 6u 2 - 2u yy + (5 2 (213) 
f 3 (y,t) = e t3 = -2[u^ - lOuuP - 5(u y ) 2 + 10u 3 ] + /3 3 



3 Construction of one-gap densities of polar op- 
erator 

In order to construct a family of one-gap densities of the polar operator we start with 
smooth one-gap solution of the KdV equation of the form || 

u(y, h) = -2p(i(y + 4Ch) + u) + 2 /s G. (314) 

Here p(z) is elliptic Weierstrass function [l^] defined in a standard manner by real pa- 
rameters ej, j = 1, 2, 3 such that e$ < e% < ei, J2j e j = 0; w ! w' are correspondingly real 
and purely imaginary semi-periods of the function p(z): p(uj a \u,uj') = e a , a = 1,2,3, 
where 

U)\ = U), L02 = w + <-<A ^3 = w , 

and C G R is arbitrary constant. 

Since in what follows we are interested in the densities of the polar operator rather 
than solutions of the Harry-Dym equation, everywhere below we omit the dynamics simply 
putting t = 0. However it should be noted that eqs. (^) have been derived from HDE 
hierarchy and completely are based on dynamics. 

The spectrum of the Schrodinger equation (|i~3| ) (with potential u(y) = u(y,0)) has 
only one gap {Ex^E^) and the spectrum boundary points have form Eq = ez + 2 /%C, 
Ei = e 2 + 2 / 3 C, E 2 = ei + 2 / 3 C. 
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To obtain all 1-gap densities of the polar operator generated by the potential ([31 

ELD 5 
R(y) 



one has to insert (314) and ( |212| ) into (^). This yields 

C 



p(iy + uj) + ±C + lf 

Now one has to check the validity of relation u = (i? -1 / 2 ) i? 1 / 2 where u is given by 

\ ' yy 

( |314| ). It turns out that it holds not for any values of C and f3\. Namely we come to the 
following 

Lemma 1. Potential u from (314) gives rise to density of the polar operator only if 
C = 2 lz C is solution of the equation: 

AC 3 - g 2 C + g 3 = 

and Pi = 0. 

Since the solution of the latter equation reads C = —e a , a = 1,2,3 one sees that 
there exist 1-gap densities only with following arrangement of spectrum boundary points: 
Eo = — e a , Ei = e2 — e a , E2 = e± — e a . So by all means one of the boundary points 
coincides with 0. 

Expressions for finite-gap densities (emerging if we put ti = in solution of the Harry 
Dym equation) read: 

r a (x) = A a [e a - p(iy + o>i + u a )], 

x = e a A a y - iA a [C(iy + u + uj a ) - r? - rj a ], a = 1,2, 3. (315) 



Here 



Ax = ^2 - , Hi = (e a - e ( g)(e Q - e 7 ), a = 1,2, 3, 



C is Weierstrass ("-function: — C'( z ) = p(z) and rj a = ((uj a ). 

Let us briefly overview functional properties of constructed densities. 
The density r a (x), a = 1,2,3 is periodic function with real period 

T a = e a A a 2\u>'\ — 2iA a 7/3, (316) 

it's spectrum consists of spectral bands [Eq,Ei] and [-E2,+°°]- 

In case a = 3 one has: £^0 = 5 E\ = e2 — e^, Ei = e\ — e^. Function r3(x) is even 
(i. e. r^{—x) = r%(x)), periodic, smooth and doesn't turn into 0. Therefore the string 
density arising in eq. (|i~2|) p% = 1 / r 2, possesses the same properties. It is density which 
M. G. Krein wrote about in his paper |I| and therefore we refer to function r%{x) as the 
Krein density. According to Krein's theorem there exists just one periodic even density 
of the string equation p = 1 / r 2 having one-gap spectrum which starts at point 0. It 
should be noted that the density r%(x) was first constructed in ||. It's period T3 first had 
been found in (IJ] in terms of elliptic Jacobi functions. It is this density which has been 
constructed here explicitly. 

In case a = 1 r\(x) is a discontinuous function. Therefore instead of smooth generating 
solution of the KdV equation (314) we take singular one, namely 

u(y,ti) = -2p( i (y + 4Cti))+ 2 / 3 C. 

Using the same technique we arrive to the density 

fi(x) = Ai[ei - p(iy + u)], 

x = eiAty - iAi[C{iy + u)) - 77], 
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As formula ( |316| ) for periods of densities r a (x) remains to be valid in this case we 
conclude that density f\(x) is a periodic even function with period T\ = 2e\ A\ \J\ —1iA\r( 
having singularities of the "cusp" type: 

2 / 

f\{x) ~ Const x ' 3 

in points x n = nT\, n £ Z. Therefore we refer to this density as cusp-periodic density. 
The boundary spectrum points are Eq = 63 — ei, Ex = e 2 — ex, E 2 = 0. 

This density is remarkable for it possesses degeneration property, namely if the first 
spectral zone shrinks into single point (Eq = E\ = —2>"f) and continuous spectrum fills 
R + (E 2 = 0) then f\ (x) turns into well-known soliton density of the Harry Dym equation: 

R 1 (y)=th 2 (^y), 

Ek{y) = -4=th(V37!/). 

In case a = 2 we obtain singular density similar to r\. Again changing the initial 
generating solution of the KdV equation by 

u{y, h) = -2p(i(y + 4Cti) + J) + 2 /a C 

we arrive to the density 

f 2 {x) = A 2 [e 2 - p(iy + oj)], 

x = e 2 A 2 y- iA 2 [((iy + u) - 77]. 

Here the properties are quite similar to the properties of the f\. Namely it has "cusps" 
at the points x = -y- +nT 2 , where T 2 is it's period with respect to x. Here the boundary 
spectrum points are Eq = e% — e 2 , E\ = 0, E 2 = e\ — e 2 . 

Note that when trying to degenerate f 2 while the first spectral zone shrinks into single 
point (Eq = E\ = —37) we get A 2 — > 00 and thus f 2 doesn't possess degeneracy property 
analogous to that of f 1 . 



4 Two-gap densities of the polar operator 

In this section we apply the technique described in sect. ^ to construct elliptic two-gap 
densities of the polar operator. As far as we know these densities weren't present in 
literature before. 

Consider the case when the Riemann surface ( |27| ) possesses a special sort of symmetry, 
namely let 

E = -y/^ + Sa, E a = -3e a + 3a, a = 1,2,3, £ 4 = v / 3& + 3a, (417) 

where a — arbitrary real parameter. The points of additional spectrum {/j^}| =1 necessary 
to define 2-gap potential of the Schrodinger equation unambiguously we choose (following 
p| ) to be the following: 



//1 = E 2 = -3e 2 + 3a, fi 2 = E 3 = -3e 3 + 3a, 
i. e. we put additional spectrum at the ends of the "main" spectrum. Then according to 



[12] the general formula for potential in terms of 2-dimensional B-functions is reduced to 
to form: 

u(y) = — 6p(iy + u>\ u, u) + 3a, (418) 
i. e. it is two-gap Lame potential. 
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In 2-gap case insertion of ( |41g| ), ( |212| ) and (213) into ( p6| ) yields 

K 



R{y) 



p 2 (iy + oj) + ap(iy + oS) + 6' 



(419) 



where 



Now as in 1-gap case one has to check the validity of the relation u(y) = (R 1 I 2 \ R 1 / 2 

V > yy 

with u being the generated 2-gap potential (418). The following statement is valid: 

Lemma 2. Potential u from ( [418] ) gives rise to density of the polar operator only if a 

is solution of the equation: 



.2 & 



^)(4a 3 -g 2 a-g 3 )=0, 



and the constants f3\ and 02 are related as follows: + fh. = 0. 

One sees that parameter a may take only 5 values: 



, /g2 

a = ±Wy,ei,e2,e3. 

This means that we get 5 two-gap densities of the form (419). Due to eqs. ( (417| ) one of 
the spectrum boundary points by all means coincides with zero. 

Each density r a (x), a = ±, 1, 2, 3 as in 1-gap cases is defined by pair of equations 



r a (x) = R a (y), x = y - / (1 - R a (y') dy 



In the case a 2 = ±a/3 (a = ±) we obtain two densities 



r±(x) 



i/3 



2 ' 



x = h -5 

a 13 - g 3 a d - g 3 



{((iy + u> - 7±) + C(% + w + 7±) - 2r]}. 



(420) 
(421) 



where 7± is a solution of the equation = — p(7±)- 

The periods of functions r±(x) implicitly described by (|420| ) and ([42 1|) are 



rp 2a , , /g7 

T± = ^ (ow \ + 2irj), a = ±J — . 

a 6 - g 3 V 3 

In the case a = e a (a = 1, 2, 3) we obtain densities 



15 . 



r a (x) 



21 



g2 



(p(iy + u) - ep){p(iy + u) - e 7 ) ' 



-y[A 



e/3 e 7 



K 



(422) 



e/3 - e 7 
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Here it is assumed that symbols {a, f3, 7} represent some transposition of numbers 
{1,2,3} 

Each constructed density r a (x), a = 1,2, 3 is a periodic function with period: 
T. _ 2A |«/| (ix _ _ 2Ai ,f (± _ _Lj , „ _ lf 2 , 3 

Among these five densities the most interesting is r + (x). It is smooth, even, nonva- 
nishing function. All of it's spectral bands [Eq,Ei], [E^-E^], [^4,+co] lie on R + and 
Eq = 0. Thus the density r + [x) is two-gap density of Krein's type Q. 

Another four densities are discontinuous functions. Their analysis is out of scope of 
this paper. 

5 Auto-Backlund transformations of finite-gap den- 
sities 

In this section we analyze the densities which can be obtained by auto-Backlund trans- 
formation of the one- and two-gap densities constructed above. Let's return for a while 



to a KdV and HD equations. Remind the following statements [11], [jj. The pair of KdV 
solutions u(y, t) and u(y, t) related by the auto-Backlund transformation can be written 
in terms of the HDE solution r(x,t) = R(y,t) as 



u = u + d y InR, 



where R corresponds to u under the relation u = [R 1 ^ 2 ) R 1 ! 2 . 

^ ' yy 

Let by f(x,t) = R(y,t) denote the HDE solution related to u in the same manner, 
namely: 

u=(R- 1 / 2 ) R 1 / 2 . (523) 
' yy 



Then as follows from [11] 



fl = |, (524) 



y 



x = y- (1-R(y',t))dy'. 







In case of decreasing density (r — > 1 when |x| — > 00) the constant b = 1 [11]. In case 
of periodic density this constant can not be obtained from the general relations and some 
additional considerations are needed. It turns out that constant b values are different for 
individual densities. Below we construct explicitly auto-Backlund transformation for 1- 
and 2-gap densities presented above and compute the value of constant b in each case. 

1-gap case 

Here we compute auto-Backlund transformation for smooth 1-gap Schrodinger poten- 
tials, u a (y) = —2p(iy + to) — e a which has been used in section || for generating the 
densities R a (y) = r a (x) given by (|315j ). 

From (HD, (HD it follows that 

u a {y) = (R^iyj) R^ h {y). (525) 
v ' yy 

The result of calculations is following. Potentials u a (y) are: 

u a (y) = -2p(iy +w + w tt ) - e a . 
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When a = 3 we obtain smooth potential which differs from its auto-Backlund trans- 
formation us(y) by shift in y. When a = 1, 2 we obtain singular potentials. Thus in 1-gap 
case the auto-Backlund transformation for the KdV equation doesn't lead out the class 
of known 1-gap potentials. 

In order to calculate auto-Backlund transformation f a (x) = R a (y) of 1-gap densities 
r a (x) = R a {y) the whole procedure of densities R a construction should be repeated 
starting from potentials u a rather than u a . The properties of densities f\{x) and ^(x) 
have been already considered in section ||. The density r^(x) differs from the smooth 
density r^(x) which possesses no roots only by shift in x by Ts /2- 

Thus we have got 5 distinct 1-gap densities of the polar operator. Densities f\{x) and 
r2{x) are continuous cusp-periodic densities. Their auto-Backlund transformations r\{x) 
r%(x) correspondingly are singular functions. Densities r$(x) and r$(x) differ only in shift 
by half-period Ti in x and represent unique smooth 1-gap Krein density whose spectrum 
starts in and spectral zones lie at R + . 

Concluding analysis of auto-Backlund transformation in 1-gap case we put expressions 
for the constant b from ( |525| ), Expressions are: 

t t ( 3 /2 e a ) 2 

b = b a = - r, a = 1,2, 6. 

(e a - e/3)(e a - e 7 J 

2-gap case 

The scheme of calculating auto-Backlund transformation of two-gap densities r a (x) = 
R a (y), a = ±,1,2,3 (see eqs. (|42Cf ) , (|422j )) is the following: one computes u a (y) by 
means of relation ( |525| ) and then use this function as generating Schrodinger potential 
in the scheme of constructing polar operator densities presented in section |3[ The result 
reads: 



r±(x) = K±(p{iy + u) ± - \Jg 2 /^ 



6 " 6 

where 



h ^ _ + =p iaK± [Q{iy + u))-rj\ 



K± 



75a 2 - 7g 2 




12a(a 3 -g 3 )a 2 ' 
Another three densities have the form: 

f a (x) = K a (p{iy + u))- ep)(p(iy + e 7 ), a = 1, 2, 3. 

x = K± + epe^j y - p'(iy + uj) - ie Q K±[((iy + u) - rj\ 



where 

K _ (75a 2 - 7g 2 ) _ _ 

2(3a2-g 2 )(12a2-g2)(fe a -^g 2 ) 

Before describing properties of the densities obtained let us note the following. Since 
we don't prove the fact that auto-Backlund transformation in periodic case preserves the 
structure of spectrum of the corresponding linear problem for the polar operator we can 
not state that spectra of the densities f±(x), r a (x) consist of only two bands. Nevertheless 
properties of these densities outlined below suggest the validity of the said hypothesis. 

Now we briefly describe properties of the obtained densities. The density f + (x) is the 
smooth periodic symmetric nonvanishing density of the Krein type. The density f_(x) 
has cusp-type singularities 



4 / 

.(x) — Const (x — Xn) <°. 
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The periodic densities f a (x) also have cusp singularities although of another type: 

2 I 

r<x{x) ? — r Const (x — x a ) /s . 

Let us note that to implement the whole scheme which led to the above expressions 
the constants f3\ and /3 2 in (|212j ) and (|213| ) should satisfy condition 



(30a/3i + /3 2 ) = 180a 2 . 

Finally calculations for the constants b in ( |525| ) for each pair of densities linked by the 
auto-Backlund transformation give: 

75a 2 - 7g 2 
* 12a(a3-g 3 )a 2 ' ° V3 
, (75a 2 -7g 2 ) 

o a = 7 — , a = e a , a = 1, 2, 3. 

2(3a 2 -g 2 )(12a 2 -g 2 )(fe Q -^g 2 
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